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Abstract. Wc perform various changes of measure in the lookdown particle system of Don- 
nelly and Kurtz. The first example is a product type h-transform related to conditioning a 
Generalized Fleming Viot process without mutation on coexistence of some genetic types in 
remote time. We give a pathwise construction of this h-transform by just "forgetting" some 
reproduction events in the lookdown particle system. We also provide an intertwining rela- 
tionship for the Wright Fisher diffusion and explicit the associated pathwise decomposition. 
The second example, called the linear or additive h-transform, concerns a wider class of mea- 
sure valued processes with spatial motion. Applications include: -a simple description of the 
additive h-transform of the Generalized Fleming Viot process, which confirms a suggestion 
of Overbcck for the usual Fleming Viot process -an immortal particle representation for the 
additive h-transform of the Dawson Watanabe process. 



1. Introduction 

Measure valued processes are usually defined as rescaled limit of particle systems. At the 
limit, the particle picture is lost. It is nevertheless often useful to keep track of the particles in 
the limiting process. First attempts to do that were concerned with a single particle, the most 
persistent one: this generated the so called spinal decompositions of superprocesses, see Roelly 
and Rouault [16], Evans [5] and Overbeck [13]. Second attempts deal with many particles, 
still the most persistent, as the infinite lineages of a supercritical superprocess, see Evans and 
O'Connell [6]. Most interesting is to keep track of all the particles; this can be achieved by the 
following trick: ordering the particles by persistence (and giving them a label called the "level" 
accordingly) allows one to keep a particle representation of the full system after taking the 
limit; the measure valued process is then represented by (a multiple of) the de Finetti measure 
of the exchangeable sequence formed by the types of the particles. This program was realized 
by Donnelly and Kurtz [5] with the construction of the lookdown particle system. 

Our aim in this article is to explain that some transformations of the law of measure val- 
ued processes, which belong to the class of /i-transforms, admit a simple interpretation when 
considered from the lookdown particle system point of view. 

We recall Doob /i-transform refers to the following operation: given a transition kernel 
Pt{x,dy) of a Markov process and a positive space time harmonic function H{t,y) for this 
kernel, meaning that: 

^ H{t,y) pt{x,dy) = H{0,x) 
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for every x and t, a new transition kernel is defined by pt{x,dy)H{t,y)/H{0,x), and the asso- 
ciated Markov process is called an /i-transform. Working with measure valued processes, we 
may choose H{t, y) to be a linear functional of the measure in which case the /i-transform is 
called additive. 

Our contribution is the following: we observe that the Radon Nikodym derivative H{t, y) / H{0, x) 
may be simply interpreted in term of the lookdown particle system in the following two cases: 

• For a probability measure valued process on {1 . . . K'} called the Generalized Fleming 
Viot process without mutation, the choice H{t, y) = e*"-^* Y[f=i for 1 < K < 
K', y a probability measure on {1, . . . , K'} and a non negative constant chosen so 
that H is harmonic, amounts to allocate the K first types to the K first particles. 
The corresponding /i-transform is the process conditioned on coexistence of the K first 
types in remote time, and the associated lookdown particle system is obtained by just 
"forgetting" some reproduction events in the original particle system. This way, we give 
a genealogy to the Wright Fisher diffusion conditioned on coexistence. We also take 
the opportunity to present an intertwining relationship for the Wright Fisher diffusion 
and explicit the associated pathwise decomposition. This adds another decomposition 
to the striking one of Swart, see [21] . 

• For a more general measure valued process on a Polish space E (incorporating mutation 
and non constant population size), the choice H{t,y) = J y{du)h{t,u) for a suitable 
function h(t, u) : M"^ x E ^ M"*" of the underlying mutation process and y a finite 
measure on E, amounts to force the first level particle to move like an /i-transform of 
the underlying spatial motion (or mutation process), and to bias the total mass process. 
This confirms a suggestion of Overbeck about the additive /i-transform of Fleming Viot 
processes, see [T^ p. 183. This also relates in the branching setting to decompositions 
of the additive /i-transforms of superprocesses found by the same author [13] using Palm 
measures. 

Our two examples, although similar, are independent: the first one may not be reduced to the 
second one, and vice versa. We stress on the change of filtration technique, learnt in Hardy 
and Harris [S], which allows us to give simple proofs of the main results. 

We first recall in Section 2.1 the lookdown construction of [3] in the case of the Generalized 
Fleming Viot process without mutation in finite state space. We look in Section 2.2 at the 
aforementioned product-type /i-transform, and prove in Section 2.3 that it may be interpreted 
as the process conditioned on coexistence of some genetic types. In Section 2.4, we compute 
the generator of the conditioned process in case the finite state space is composed of only two 
types, and recognize it as the generator of a Generalized Fleming Viot process with immigration. 
This Section also contains the statement and the interpretation of the intertwining relationship. 
Section 3.1 starts with the introduction of a lookdown construction allowing for mutation and 
non constant population (also extracted from [3]). We then present in Section 3.2 the additive 
/i-transform of the associated measure valued process. Section 3.3 is concerned with applications 
in two classical cases: Dawson Watanabe processes and Generalized Fleming Viot processes. 

2. A PRODUCT TYPE /l-TRANSFORM 

2.1. The construction of the Generalized Fleming Viot Process without mutation. 

Donnelly and Kurtz introduced in [3] a population model evolving in continuous time. We 
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present in this Section the particular case of the Generahzed Fleming Viot (GFV) process 
without mutation in which we are interested. A more general framework will be introduced in 
Section 3. 

We denote by V^o the space of partitions of the set of integers N = {1, 2, 3, . . .}. We assume 
c > 0. We define jJ^ the measure on Voo assigning mass one to partitions with a unique non 
trivial block consisting of two different integers, and call it the Kingman measure. We let v be 
a measure on (0, 1] satisfying J^^ x^fi^dx) < oo. We denote by dt the Lebesgue measure on 
M+, and by px the law of the exchangeable partition of N with a unique non trivial block with 
asymptotic frequency x: If (f/j)^^^ is a sequence of independent Bernoulli random variables with 
parameter x, then the partition vr whose unique non trivial block contains the integers i such 
that Ui = 1 has law p^. Finally, we define N{dt,dn) the Poisson point measure on x Voo 
with intensity 



Let Rq be a random probability measure on the finite state space E = {1,2, . . . , K'} for 
K' > 2. Assume Rq is independent of A^. Conditionally on {Ro,N), we define the lookdown 
particle system X = {Xt{n),t > 0,ra G N), and we call Xt{n) the type of the particle at level 
n e N at time t > 0: 

• The initial state (Xo(n),n G N) is an exchangeable sequence valued in E with de 
Finetti's measure Rq: That means that, conditionally on Rq, (Xo(n),n G N) form a 
sequence of independent random variables with law Rq. 

• At each atom {t, vr) of A^, we associate a reproduction event as follows: let ji < j2 < • • • 
be the elements of the unique block of the partition vr which is not a singleton (either 
it is a doubleton or an infinite set). The individuals ji < j2 < • • • at time t are declared 
to be the children of the individual ji at time t—, and receive the type of the parent ji, 
whereas the types of all the other individuals are shifted upwards accordingly, keeping 
the order they had before the birth event: for each integer i, Xt{j() = Xt_{ji) and for 
each k ^ G N}, Xt{k) = Xt-{k - #4) with := {£ > < k} and #4 the 
cardinal of the set J^. 

• For each G N, the type Xt{n) of the particle at level n do not evolve between repro- 
duction events which affect level n. 

Remark 2.1. The integrability condition J^^ < oo ensures that finitely many repro- 

duction events change the type of a particle at a given level in a finite time interval, ensuring 
the preceeding definition makes sense. 

For a fixed t > 0, the sequence {Xt{n),n G N) is exchangeable according to Proposition 3.1 
of [3] . This allows to define the random probability measure Rt on E as the de Finetti measure 
of the sequence {Xt {n),n G N) : 



and this defines a random process indexed by the set of non negative real numbers M"*" since 
a cadlag version of the process (where the space M.f{E) of finite measures on E is endowed 
with the topology of weak convergence) is shown to exist in [3]. The process R is called 
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the Generahzed Fleming Viot process (GFV process for short) without mutation. We stress 
that, conditionally given Rt, the random variables {Xt{n),n & N) on E are independent and 
identically distributed according to the probability measure Rt thanks to de Finetti's Theorem. 
We will denote by P the law of X. We now introduce the relevant filtrations we will work with: 

• (J-f = (r{{Xs{n),n G N), < s < t}) corresponds to the filtration of the particle system. 

• (Gt = (^{Rs, < s < t}) corresponds to the filtration of the measure-valued process R. 

Notice that X is a Markov process with respect to the filtration J-", and that i? is a Markov 
process with respect to the filtration Q. 

2.2. A pathwise construction of an /i-transform. 

2.2.1. Results. The proofs of the results contained in this Subsection may be found in 12.2.21 
Fix 1 < K < K' . We assume from now on and until the end of Section |2] that: 

K 

(2) E(J]i?o{2})>0, 

i=l 

to avoid empty definitions in the following. Recall the definition of the particle system X 
associated with R in Section 12.11 We define from X a new particle system X^ as follows: 

(i) The finite sequence (Xg (j), 1 < j < is a uniform permutation of {1, . . . , K}, and, 
independently, the sequence (Xo(j), j > K + 1) is exchangeable with asymptotic fre- 
quencies Rq , where R^ is the random probabihty measure with law: 

(ii) The reproduction events are given by the restriction of the Poisson point measure N 
(defined as in Subsection 12.11) to V := {(s, vr), 7r|[i^] = {{1}, {2}, . . . , {f^}}}, where t^\[k] 
is the restriction of the partition vr of N to {1, . . . , i^}, that is the atoms of N for which 
the reproductions events do not involve more than one of the first K levels. 

Remark 12.11 ensures that this definition of the particle system X^ makes sense. 

Remark 2.2. Note that the particle system (XQ(j), j > l) is no more exchangeable due to the 
constraint on the K initial levels. Nevertheless, the particle system (Xq(j),j > is still 
exchangeable, and we shall view the K first levels as K independent sources of immigration in 
forthcoming Section [2^ 

We also need the definition of the lowest level L{t) at which the first K types appear: 

(3) Lit) = inf{z > if, {1, . . . , i^} C {Xi(l), . . . , Xi(z)}}, 

with the convention that inf{0} = oo. The random variable L{{]) is finite if and only if 
nj=i -Ro{0 > O5 P-a.s., thanks to de Finetti's Theorem. The process L{t) is Ft measurable, 
but not Qt measurable. Notice the random variable L{t) is an instance of the famous coupon 
collector problem, based here on a random probability measure. We define, for z > 1, the 
pushing rates Vi at level i: 

r, = '^'~^K + [ u{dx) {1 - {1 - xY - ix{l - xy~^) . 
2 i(o,i] 
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Notice that ri = and that rj is finite for every i > 1 since J^^ v{dx) < oo. From the 
construction of the lookdown particle system, these pushing rates may be understood as 
the rate at which a type at level i is pushed up to higher levels (not necessarily 2 + 1) by 
reproduction events at lower levels. Let us define a process Q = {Qt,t > 0) as follows: 

^* P(L(0) =K) ■ 
Lemma 2.3. The process Q = {Qt,t > 0) is a non negative T -martingale, and 

(4) VA G Tt, nx'' eA)=E (1a(X) Qt) . 

By projection on the smaller filtration Qt, we deduce Lemma 12.41 We need the following 
definition of the process: 

E(nf=i ^o{^}) 

Lemma 2.4. The process M = {Mt,t > 0) is a non negative Q-martingale. 

This fact allows to define the process = {R^,t > 0) absolutely continuous with respect 
to R = {Rt,t > 0) on each Qt, t>0, with Radon Nykodim derivative: 

(5) VA G Qt, P(i?^ eA)=E {1a{R) Mt) . 
The process R^ is the product type /i-transform of interest. 
Remark 2.5. Definition ([5]) agrees with the definition of Rq . 

Remark 2.6. Intuitively, the ponderation by M favours the paths in which the K first types 
are present in equal proportion. 

We shall deduce from Lemma 12.31 and Lemma 12.41 the following Theorem. 

Theorem 2.7 . Let I < K < K' . We have that: 

(a) The limit of the empirical measure: 

1 ^ 

R^idx) := lim — > 5Yh(„\{dx) 

n=l 

exists a.s. 

(b) The process {R'l,t > 0) is distributed as {Rf ,t > 0). 

Let us comment on these results. The process X'^ is constructed by changing the initial 
condition and forgetting (as soon as K > 2) specific reproduction events in the lookdown 
particle system of X. Lemma 12.31 tells us that this procedure selects the configurations of 
X in which the K lowest levels are filled with the K first types at initial time without any 
"interaction" between these K first levels at a further time. Theorem 12 . 71 tells us that the process 
R^ constructed in this way is an /^-transform of R and Lemma 12.41 yields the following simple 
probabilistic interpretation of the Radon Nikodym derivative in equation ([S]): the numerator 
is proportional to the probability that the K lowest levels are occupied by the K first types at 
time t, whereas the denominator is proportional to the probability that the K lowest levels are 
occupied by the K first types at time 0. We shall see in Section [2l3] that the processes X^ and 
R'' also arise by conditioning the processes X and R on coexistence of the first K types. 



6 OLIVIER HENARD 

2.2.2. Proofs. 

Proof of Lemma \2.'d From the de Finetti theorem, conditionally on i?^, the random variables 
{Xt{i),i G N) are independent and identically distributed according to Rt- This implies that: 

K 

(6) F{L{t) = KlGt) = K\ l[Rt{i}. 

1=1 

In particular, we have: 

K 

F{L{0) = K) = K\ E([[Ro{i}), 

1=1 

which, together with ([2]), ensures that Qt is well defined. 

Then, let us define W = {tt, tt\ik] = {{!}, {2}, . . . , {K}}}, and 1^ = {(s, vr), < s < t, vr G 
W}, and also the set difference W" = Voo\W and V^^ = {(s,7r),0 < s < t,7i e W}. We 
observe that: 

• from the de Finetti Theorem, the law of Xq, as defined in (i), is that of Xq conditioned 
on {L(0) = K}. 

• The law of the restriction of a Poisson point measure on a given subset is that of a Poisson 
point measure conditioned on having no atoms outside this subset: thus conditioned 
on having no atoms in (this event has positive probability) is the restriction of to 
Vt. 

Since the two conditionings are independent, we have, for A G J^t'- 

¥{X^ eA)= F{X G A\{L{0) = K}n {N{Vt') = 0}) 

l{L(o)=x}n{Af(yi'=)=o} 



^ (^l^(^) p(^(o) ^ K)¥{N{Vf) = 0) 
We compute: 

^{W") = fi^W) + I u{dx)p^{W 

K{K-1 



'(0 

c + [ i^{dx) {1 - {1 - x f - Kx{l - x)^-^) 

i(0,l] 



2 

= rx- 

This implies from the construction of that: 

(8) P(A^(^/) = 0) = e-'^^^')* = e"^^*. 
Notice that 

(9) {L(t) = K} = {L(0) = K} n {NiVn = 0}. 
From (IT]), dH]) and ([9]), we deduce that: 



□ 
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Proof of Lemma 2. 4' We know from Lemma 12.31 that {Qt,t > 0) is a J-"- martingale. Since 
Gt C J^t for every t > 0, we deduce that {E^QtlGt),^ > 0) is a ^-martingale. Now, we notice 
that: 

using ([6]) for the second equality. 

□ 

Proof of Theorem 2/l_. From Lemma [2.31 is absolutely continuous with respect to X on J^f 
The existence of the almost sure limit of the empirical measure claimed in point (a) follows 
from ([T]). We now project on Gt the absolute continuity relationship on J^t given in Lemma [231 
Let A e Gt- 

P(i?'^ G A) = E {lA{R)Qt) = E {lA{RmQt\Gt)) = E {lA{R)Mt) = e A), 

where we use Lemma 12.31 for the first equality and the definition of R^ for the last equality. 
This proves point (b). □ 

2.3. The /i-transform as a conditioned process. We study the conditionings associated 
with the /i-transforms. 

Let 1 < K < K' . Assumption ([2]) allows us to define a family of processes on F by: 

K 

MA e Gt. HR^-'^ e A) = g A| ^ o), 

i=l 

and the associated particle system X*^-*^ on G by: 

K 

\/A G J^t, P(X(^*) eA) = ¥{X G A| ^ 0). 

i=l 

The process i?*^-*^ thus corresponds to the process R conditioned on coexistence of each of the 
first K types at time t. It is not easy to derive the probabilistic structure of the particle system 
^(>*) all J^f Nevertheless, for fixed s > 0, the probabilistic structure of X*^-*^ on the sigma 
algebra J^s simplifies as t goes to infinity. 

The following Theorem may be seen as a generalization of Theorem 3.7.1.1 of Lambert [TT] 
which, building on the work of Kimura [S], deals with the case = 0. We write for the law 
of L (defined in ([3])) conditionally on {^(0) = i}. For / an interval of M"^ and F a process 
indexed by M"*", we denote by Fi the restriction of F on the interval I. 

Theorem 2.8. Let s > be fixed. Assume that 

(10) PWLW < PC) ^ 

^ ^ t^oo ^kW) < OO) 

We have: 

(i) The family of processes X^^^^ weakly converges as t —t- oo towards the process . 

(ii) The family of processes *j weakly converges as t — t- oo towards the process . 

Lemma [2.131 gives a sufficient condition for ([TOl) to be satisfied. 
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Remark 2.9. The case K = 1 corresponds to a non degenerate conditioning since the event 
{-R({1} 7^ for every t} has positive probabihty under 

Remark 2.10. Assume K > 2. The following property 

(GDI) P(inf {t > 0, L{t) = cx)} < cx)) = 1, 

is independent of K (used to define L in ([3])). (GDI) property is easily seen to correspond to 
the Goming Down from Infinity property for the A-coalescent associated with the GFV process 
R (whence the acronym (GDI)). We refer to Schweinsberg [19] for more details about this 
property. We conjecture (GDI) property is equivalent to our assumption ( ITU]) . We were unable 
to prove it. 

Remark 2.11. It should still be possible to interpret the processes X'^ and R^ as conditioned 
processes, without assuming ffTOj) . In that more general case, should correspond to X 
conditioned by the event {lim inft_^oo J^^^^ i?t{z} > 0} (which has null probability as soon as 
K > 2). 

Proof. First observation is that, from the Kingman's paintbox construction for exchangeable 
random partition, we have: Y[f=i -^tl^} 7^ if and only if L(t) < oo, P a.s. This gives, for any 
A e J",: 

K 



r(An{Ut,Rt{^}^o} 

FiA\l[Rt{i}^0^ - ^ 



. " nnL^tW^o) 

_ ¥{An{L{t) < oo}) 
P(L(t) < oo) 

Now, using the Markov property, we have: 

¥{A n {L{t) < oo}) = ¥{A n {L{s) = K} f] {L{t) < oo}) + ¥{A n {L{s) > K + 1} n {L(t) < oo}) 
= E(lAn{L(s)=x} - s) < oo)) + E(l^n{L(.)>i^+i} '^L{s){.L{t - s) < oo)). 



where L is an independent copy of L. 

Let £ G N. We can couple the processes L under P^ and L under P^+i on the same lookdown 
graph by using the same reproduction events. More precisely, imagine that we distinguish the 
particles at level I and £ + 1 at initial time, giving each of them a special type shared by no 
other particles. Then the first levels at which these two types may be found at time t yields 
a coupling of L{t) under P^ and L{t) under P^+i. Let us denote by (L£,L^+i) this coupling: 
is distributed as L under P^ and L^+i is distributed as L under P^+i. By the ordering by 
persistence property of the lookdown graph, we have that, for every t > 0: 

U{t) < L,+i(t), 

whence: 



(11) 



P,+i(L(t) < oo) <P,(L(t) < oo) 
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for every integer l. Therefore, we have: 

nUrML{s)>K+i}^L{s){L{t -s)< cx))) < P(A n {L{s) >K + l]) ^k+iW - s) < cx)) 
Our assumption ( ITU]) now imphes: 

P(An < CX)}) 



P/^(L(t - s) < oo) t->oo 
Setting A = Vt, this also yields: 

P(L(t) < oo) 



^ P(An {L(s) = K}). 



^ X X niis) = K). 

Taking the ratio, we find that: 

P(An{L(t) < oo}) P(A n {L(s) = K}) 

P(L(t) < oo) t^oo P(L(s) = K) 

We also have that P(L(s) = ii") = P(L(0) = K) e~^'^^ since Q is a ^-martingale from Lemma 
12.31 Altogether, we find that: 

lim P(A| ni?.{^} ^ 0}) = E (lA(^)p^l^e-^) = P(X'^ G A) 

where the last equality corresponds to Lemma [2731 This implies the convergence in law of X(^*) 
towards as t — ?■ oo. This proves the first point. 

The proof of (ii) is similar to the one for (i). □ 

Remark 2.12. Having introduced in the preceeding proof the coupling [Lk, Lk+i), we may 
complete the Remark 12. lOt It is possible to prove that, if (GDI) holds and for each t >0, 

(j — F{LK+iit) < oo|Lx(t) < j)) is non increasing 

then din]) holds. 

We now give a sufficient condition for f fTO]) to be satisfied. 
Lemma 2.13. IfYjj>KV' ^ ^^^^ ( fiOj) holds. 
Proof. A lower bound for Pi^(L(t) < oo) is easily found: 

(12) e""^* = fK{L{t) = K)< PA'(L(t) < oo). 

We now look for an upper bound for FK+i{L{t) < oo). Recall the non decreasing pure jump 
process L jumps with intensity rj when L = j. 
We may write, under Fk+i'- 

sup {t, L{t) < oo} = Tj 

j>K+l 

where, conditionally given the range of the function L we shall denote by L^^^ < L^^"^ < . . ., 
the sequence {Tj,j > A' + 1) is a sequence of independent exponential random variable with 
parameter r^. Since ('"j)j>a'+i forms an increasing sequence and the function L has jumps 
greater or equal to one, we have for each j > K + 1, 

(13) r^, > Tj. 
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Let [Tj , j > i^' + 1 ) be a sequence of independent exponential random variables with parameter 
('^j, j > K). We compute, for < A < r^+i '■ 

j>K+l 

<nYl > ^) 

j>K+l 

= P(exp (A Tj) > exp (At)) 

j>K+l 

< exp (-At) E(exp (A ^ Tj)) 

j>K+l 

= exp(-A*) n 7^ 

j>K+l ' 1 ^ 

= exp(-At+ log(l + ;r4A)) 

j>K+l ^ 

< exp 



— A ' 

j>K+l 



where we use f[T5]) for the first inequality and the Markov inequality for the second inequality. 
From the assumption, X]j>ii'+i V''^i is finite, which implies also that X]j>A'+i finite. 
Taking A = {rx + rK+i)/2, we obtain that: 

(14) fK+iW) < oo) < Cexp (^-!iL±I^t 



for the finite constant C = exp A Ylj>K -'-/(^i ~ associated with this choice of A. Altogether, 
using (|T2|) and f|T^ . we have that: 

- Pi^(L(t) < oo) - ^ V 2 
Letting t tend to oo, we get the required limit. □ 



As an immediate corollary, we get the following result, which ensures that ( iTOl) is satisfied in 
the most interesting cases. 

Corollary 2.14. // (c 7^ 0), or (c = and u[dx) = x~^~°(l — x)""^l(o,i)(x)(ix for some 
1 < a <2), then (E^j holds. 

Proof. If (c 7^ 0), > cj{j — l)/2, and thus ^j>xl/Tj < 00. If (c = and 
i^((ia;) = x"^~°(l — x)°~^ 1(0,1) (x)(ix for some 1 < a < 2), using the equality: 

Tj+i — r j = / j(l — xy~^x^v{dx)] 

J {0,1] 
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one may find at Lemma 2 of Limic and Sturm [12] , we get 

r,+i-r, =jBeta(2-a,j + a-l) ~ r{2 - a)r~\ 

wfiere, as usual, Uj ~ Vj means that Vj ^ for j large enough and Xmxj^oaUjjvj = 1. We 
conclude that rj ~ r(2 — a)j"/a, and then: J2j>K -'-Aj < since 1 < a < 2. 

j—^oo -' — 

Lemma 12.131 allows to conclude that (fTOj) holds in both cases. 

□ 

2.4. The immigration interpretation. We develop further the two following examples: 

(i) K = K' = 2: this amounts (provided condition ffTOl) is satisfied) on conditioning a 
two-type GFV process on coexistence of each type. 

(ii) 1 = K < K' = 2: this amounts (provided ffTOj) is satisfied) on conditioning a two-type 
GFV process on absorbtion by the first type. 

We regard the K[= 1 or 2) first level particles in X'^ as K external sources of immigration in 
the population now assimilated to the particle system {X^{n),n > K + 1) and decompose the 
generator of the process R'^ accordingly. We refer to Foucart [7] for a study of GFV processes 
with one source of immigration [K = 1 here). 

Since K' = 2, the resulting probability measure-valued process R = {Rt,t > 0) and R^ = 
{Rfy t > 0) on {1, 2} may be simply described by the [0, l]-valued processes -R{1} = {Rt{l}, t > 
0) and R^{1} = (i?f{l},t > 0) respectively. For the sake of simplicity, we will just write R for 
R{1} and R'^ for i?''{l} respectively. 

We recall that the infinitesimal generator of R is given by: 

Gf{x) = lcx{l-x)f"{x)+x [ u{dy)[f{x{l-y)+y)-f{x)] + {l-x) [ u{dy)[f{x{l-y))-f{x)] 

^ ^(0,1] -^(0,1] 

for all / G C^([0, 1]), the space of twice differentiable functions with continuous derivatives, and 
X G [0, 1], see Bertoin and Le Gall [1]. 

2.4.1. We assume K = K' = 2. We define, for / e C2([0, 1]), and x G [0, 1]: 

G'fix) = c(l - 2x)f{x) + / y{l- yMdy)[f{x{l - y) + y) - f{x)] 

J{0,1] 

+ [ y(l-y)i,(dy)[f{x{l-y))-f{x)l 

J{0,1] 

and 

G'fix) = lcx{l - x)f{x) + x [ (1 - y)Mdy)[f{x{l - y) + y)) - /(x)] 
^ i(o,i] 

+ (l-x) / (l-y)V(rfy)[/(x(l-y))-/(x)]. 

7(0,1] 

Proposition 2.15. Assume K = K' = 2. The operator + G^ is a generator for R^\ 
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Remark 2.16. When the measure u is null, the process R is called a Wright Fisher (WF in 
the following) diffusion. In that case, the process R^ may be seen as a WF diffusion with 
immigration, where the two first level particles induce continuous immigration (according to 
G^) of both types 1 and 2 in the original population (which evolves according to = G in 
that case). 

When the measure u is not null, the process R^ may still be recognized as a GFV process 
with immigration, but the generator G^ is no more that of the initial GFV process G: the two 
first level particles induce both continuous and discontinuous immigration (according to G'^) of 
types 1 and 2 in a population with a reduced reproduction (the measure ^'{dy) is ponderated 
by a factor (1 — y)^ < 1 in G^). 

Proof. Let us denote by G^ the generator of R^. The process R^ is the Doob /i-transform of R 
for the following function H, which is space time harmonic according to Lemma 12.41 

H{t,x) = x(l -x)e"^*. 

From the definition of the generator, for / G C^([0, 1]), and x G [0, 1]: 

f{Rt)H{t,Rt) - f{Ro)H{0,Ro) - [ ds G{H{s, .)f){Rs) - [ ds dtH{., Rs){s)f{Rs) 



is Q martingale, where in the first integrand G acts on a; — )■ f{x)H{s,x). Therefore, on 
{H{0,Ro) ^ 0}, the process 



H(0,Rti)-" " J„ //(0,B„) H(a,R,) J„ ff(0,fl„) H(a,R 

is again Q martingale under P. This implies that: 

m ) - m.) - 1 ds -^^^^ — I ds ^^^/(^J 

is a ^ martingale under P. We thus have: 

(15) GVM^(^^M±p^)w 

a formula which holds for any t > 0. A simple computation completes the proof of the propo- 
sition. □ 

Using the particle system X^, we also have the following intuitive interpretation of the 
generator G^ in the case of a pure jump GFV process (c = 0). Let us decompose the intensity 
measure u as follows: 

u{dy) = 2y{l - y)u{dy) + (1 - yfu{dy) + y^u{dy). 

(1) The first term is the sum of the two measures y{l — y)v{dy) appearing in each integrand 
of the generator G*^ and each of these measures corresponds to the intensity of the 
reproduction events involving level 1 and not level 2, or level 2 and not level 1 (these 
events have probability y{l — y) when the reproduction involves a fraction y of the 
population). 
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(2) The second term is the measure (1 — yYv{dy) appearing in the generator and cor- 
responds to the intensity of the reproduction events involving neither level 1 nor level 2 
(this event has probability (1 — yY when the reproduction involves a fraction y of the 
population). 

(3) The third term does not appear in the generators G° and G^: it corresponds to the 
intensity of the reproduction events involving both level 1 and 2, and these events have 
been discarded in the construction of X^. 

2.4.2. We assume K = 1, K' = 2. Note that the case K = 1 differs from the case K = 2, 
since the event {Rt 7^ for every t} has positive probability under ([2]). Let us define, for 
/ e C2([0,l]), and X e [0,1]: 

I^fix) = c(l - + / yu{dy)[f{x{l - y) + y) - f{x)] 

Jm 

and 

I' fix) = lcx{l - x)f"{x) +x [ (1 - y)u{dy)[f {xil - y) + y)) - f{x)] 

+ (1 -x) [ (1 - y)uidy)[f {xil - y)) - f{x)]. 

J{0,1] 

We can then prove the analog of Proposition 12.151 in that setting. 

Proposition 2.17. Assume K = 1, K' = 2. The operator + is a generator for the Markov 
process . 

In particular, we recover the well known fact that a WF diffusion conditioned on fixation at 
1 (that is, Rt = 1 for t large enough) may be viewed as a WF process with immigration, see [1] 
for instance. 

Proof. The proof is similar to that of Proposition 12.151 Here we use an /i-transform with the 
function 

H(t, x) = X. 

This function is space time harmonic according to Lemma [2.41 (recall ri = 0). □ 

Here again, we may have obtained the generator of R^ without computation. We explain the 
procedure in the case of a pure jump GFV process (c = 0). We decompose u as follows: 

iy{dy) = yu{dy) + (1 - y)i^idy). 

(1) The first term is the measure yv[dy) appearing in the generator P . This is the intensity 
of the reproduction events involving level 1 particle. We interpret them as immigration 
events. 

(2) The second term is the measure (1 — y)v{dy) appearing in the generator . This is the 
intensity of the reproduction events not involving level 1 particle. We interpret them as 
reproduction events. 

(3) Summing the two measures yv{dy) and (1 — y)v{dy), we recover this time the full 
measure z/((iy) since no reproduction events are discarded in the case K = 1. 
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2.4.3. Intertwining. We assume K' = 2 and u = (for the sake of simphcity). In this Subsec- 
tion, we observe that the processes L and R are intertwined in the sense of Rogers and Pitman 
[T7]. This adds another decomposition to the striking one of Swart, see [2T], which does not 
admit a clear interpretation from the lookdown particle system. 

In fact, we find it more convenient to prove rather that and R are intertwined, where: 

L\t) = mi{z > 1,1 G {Xt{l),...,Xt{i)}} 

is the first level occupied by a type 1 particle. The process is valued in NU {oo}, and jumps 
by 1 at rate c ^^^~^^ when at i. In particular, 1 is an absorbing point for L^. Notice also that 
the process R'^ studied in Subsection 12.4.21 is the process R conditioned on {L^ = !}• 
Let us define the following kernel: 

K{x, i) = il- x)^-ixl(o,i](x) + l|oo} Wl{o}(a:), X G [0, 1], ^ G N U {oo}. 
acting on function f{x,t) as follows: 

kf{x)= J2 K{x,l)f{x,l). 

i<e<oo 

We slightly abuse of notation by still denoting by G the generator of the Wright Fisher diffusion: 

G/(x) = ^cx(l-x)r(x) 

acting on / G C^([0, 1]). We denote G the generator defined for £ < oo by: 

Gf{x, i) = icx(l-x)9,,/(x, i)+c [(1 -x)-ii- l)x] dj{x, i)+c^-^^^ [/(x, ^ + 1) - /(x, i)] 

and for £ = oo by: Gf{x, oo) = 0. This generator acts on functions / such that /, as a function 
of X, belongs to C^([0, 1]). The intertwining relationship reads as follows. 

Proposition 2.18. Let f be in the domain of G and x G [0, 1]. The kernel K intertwins the 
generators G and G in the sense that: 

KG{f){x) = GK{f){x). 

The proof consists in a long but simple calculation and is eluded. A similar intertwining 
relation also holds for u 0, but the generator G is then more complicated (because and R 
may jump together in that case). The intertwining relation implies that the first coordinate of 
the process with generator G is an autonomous Markov process with generator G. 

Now, it may be justified that the generator G is the generator of {R, L^). In fact, the process 
is plainly Markov in its own filtration and jumps from £ to £ + 1 at rate ci{i — l)/2. Then 
conditionally on the value of L = i, we view the i first particles as £ immigrants, and the process 
i? as a Wright Fisher diffusion with i sources of immigration, i — 1 sources of type 2 and one 
source of type 1, whence the drift term c [(1 — x) — {i — l)x] thanks to similar calculations as 
in 12.4.11 We thus obtain the following pathwise decomposition of a Wright Fisher diffusion, 
which is another way to express the intertwining relation: 

• Conditionally on {Rq = x}, the initial value L^{0) has law: 

F{L\0) = = (1 - xY~'x + l{oo}(^)l{o}(a:), i > 1- 
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• Conditionally on [Rq, L^{0)), the process is a pure jump Markov process, which 
jumps from £ to £ + 1 at rate c£{i — l)/2 if £ < oo. 

• Conditionally on [Rq, L^), the process i? is a Wright Fisher diffusion with immmigration, 
with generator given by: 

^cx{l - x)f'\x) + l{ii<oo} c [(1 - x) - (L^ - l)x] fix). 

3. The additive /i-transform 

In this Section, we derive another example of an /i-transform (of measure valued processes) 
admitting a simple construction from the lookdown particle system. 

3.1. The general construction of the lookdown particle system. We first present a 
more general construction of an exchangeable particle system, which allows to deal with type 
mutation and non constant population size. We recall this model was defined (in even greater 
generality) in [3]. 

Let E he a. Polish space. We consider a triple {Ro,Y,U) constructed as follows. Rq stands 
for a probability measure on E, and Y = {Yt,t > 0) and U = {Ut,t > 0) for two non negative 
real valued processes. We assume that Uq = and U is non decreasing, so that U admits a 
unique decomposition Ut = + X]s<t '^^s where U'' is continuous (with Stieltjes measure 
denoted by dU'^) and AUs = Us — Ug-- We assume that is an absorbing point for Y , and 
set t{Y) = inf {t > 0,Yt = 0} the extinction time of Y. We also assume that for each s > 0, 
AUs < Y^- Conditionally on U and Y, we define two point measures A^'' and A^'^ on IR+ x Voo, 
where Voo denotes the set of partition of N: 

• A^'' = Yl ^{t,TT){dt,d'n') where the exchangeable partitions tt of N are indepen- 

0<t<T{Y),AUt¥=0 

dent and have a unique non trivial block with asymptotic frequency y/AUt/Yf. 

• N'' = ^ 6(^t,n){dt,d7r) is an independent Poisson point measure with intensity 

o<t<r(y) 

{dU^/lYt)"^) X fi^, and the Kingman measure fi'^ assigns mass one to partitions with 
a unique non trivial block consisting of two different integers, and mass to the others. 

Conditionally on {Rq, Y, U), we then define a particle system X = {Xt{n), <t < r(F), n G N) 
as follows: 

• The initial state {XQ{n),n G N) is an exchangeable sequence valued in E with de 
Finetti's measure -Rq- 

• At each atom (t, vr) of N := N'' + iV, we associate a reproduction event as follows: 
let ji < j2 < ■ ■ ■ he the elements of the unique block of the partition vr which is not a 
singleton (either it is a doubleton if (t, vr) is an atom of A^*^ or an infinite set if (t, vr) is 
an atom of A^^). The individuals ji < j2 < • • • at time t are declared to be the children 
of the individual ji at time t—, and receive the type of the parent ji, whereas the types 
of all the other individuals are shifted upwards accordingly, keeping the order they had 
before the birth event: for each integer i, Xt{ji) = Xt^{ji) and for each k ^ {ji, i & N}, 
Xt{k) = Xt^ik — 4^Jk) with Jfc := {i > < k} and #Jfc the cardinal of the set Jfc. 
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• Between the reproduction events, the type Xt{n) of the particle at level n mutates 
according to a Markov process with cadlag paths in E, with law (P^jX G E) when 
started at x G -E, independently for each n. 

This defines the particle system X on [0,t{Y)). Since the process Xs{j) admits a limit as 
s goes to t(Y) for each j, we may set Xt{j) = lims_^T-(y) ^^(j) for each t > t(Y) and this 
sequence is still exchangeable according to Proposition 3.1 of [B]. Conditionally on (.Ro,^,^^), 
the sequence {Xt{n), G N)) is well defined for each t G M"*" and exchangeable according to the 
same Proposition 3.1 of [3]. We denote by Rt its de Finetti measure: 

1 ^ 

Rt{dx) = lim — \2Sx,in){dx), 

n=l 

and this defines a random process indexed by the set of non negative real numbers M'^ since a 
cadlag version of the process R is shown to exist in [3] . We finally define the measure valued 
process of interest Z by: 

(16) {Zt,t>0) = (YtRt,t>0). 

The finite measure Z represents the distribution of a population distributed in a space E, the 
process Y corresponds to the total population size, and U tracks the resampling inside the 
population. We stress that, conditionally given R^, the random variables {Xt{n),n G N) on 
are independent and identically distributed according to the probability measure Rt thanks to 
the de Finetti Theorem. 

We will denote by P the law of the triple [Y, U,X). We introduce the relevant filtrations: 

• (J-f = cr((Ys, s < t), {Xg, s < t))) corresponds to the filtration of the particle system and 
the total population size. 

• {Gt = o-{Zs, s <t)) corresponds to the filtration of the resulting measure valued process. 

• Vt is the filtration induced by the canonical process under P. 

We shall use the classical notation = J fi{dx)f{x) for a non negative map f : E ^ M. and 
/i G A^/. Note that Yt = Zt{l), and thus Y is ^-measurable. 

3.2. The additive /i-transform. We denote by A^^ the space of finite measures on E. We 
call a non-negative function H on [0, oo) x A4f a. space-time harmonic function for P when the 
process {H{t, Zt),t > 0) is a martingale under P. The /i-transform Z^ of Z associated with H 
is then defined by: 

(17) ,Ae5..,nz''.^A)^:J^ll^nZeA). 

for every t > 0. Furthermore, an /i-transform is called additive if there exists a non-negative 
function {ht{x),t > 0,x & E) such that H(t, Zt) = Ztiht). 

Remark 3.1. Loosely speaking, an additive /i-transform of the form ( ITTl) favours the paths for 
which the population (represented by the measure valued process) is large where h is large. 
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3.2.1. Statement of the results. Let ^ be the canonical process under P^-. We assume that 
{Yt/m{t),t > 0) and {m{t)ht{^t),t > 0) are martingales in their own filtrations for some positive 
deterministic function m. We assume from now on that 

E{YoRo{ho)) > 0. 
Under this assumption, we define (the law of) a new process 

(F^[/^x'^) 

by the following requirements: 

(i) The initial condition satisfies: 

(ii) Conditionally on {Yq,Rq), and provided i?o(/io) > 0, -^o(-'-) distributed according to: 

WA e I?o,P(Xo^l) e A\R'o = /i) = E (^MgMi^(Xo(l))|i?o = fi 

and (Xq(?7,),?7, > 2) is an exchangeable random sequence with de Finetti's measure -Rq. 

(iii) Conditionally on (Yq,Rq,Xq{1)), the process {Y^,U^) is distributed according to: 

(18) e Gt, P((F^ U^) e A\Y^ = x) = E (YlV^ u)\Y, = x 



X m(t) 

(iv) Conditionally on {Y^,U^ ,Rq,Xq{\)), X^{\) is distributed according to: 
(19) \/A G A, nx\l) e A\X^il) = x)=E fMMll !!^i^(x(l))|Xo(l) = X 



h(j{x) m(0) 

(v) The rest of the definition of X^' is the same as the one given for X, namely: 

— for n >2, between the reproduction events, the type Xl^{n) of the particle at level 
n mutates according to a Markov process in E with law (P^, x E E) when started 
at x E E, independently for each n. 

— at each atom (t, tt) oi N = N'' + N^, with N'' and derived from and , a 
reproduction event is associated as previously. 

Note that the law of the initial condition Zq specified by (i) is different from that of Zq only 
for random Zq. Also, notice that items (iii) and (iv) are meaningful since both (Yt/m{t),t > 0) 
and {m{t)ht{Xt{l)),t > 0) are assumed to be martingales. Last, we observe from f|T8l) that 
P(y/' = 0) = for each t > 0, which implies P(r(F'') = oo) = 1 since is assumed to be 
absorbing. We will assume that {Y, U, X) and [Y'^, U^, X^) are defined on a common probability 
space with probability measure P, and denote the expectation by E. 
Let us define a process S = {St, t > 0) by: 

_ htiXtjl)) Yt 

Lemma 3.2. The process {S = St,t > 0) is a non negative T -martingale, and 
(20) VA G J-i, P(X^ G A) = E {1a{X) St) . 
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We then define the process T: 



E{Zo{ho)) 

Using Lemma [3. 2 [ and projecting on the filtration Qt, we deduce Lemma [3.31 

Lemma 3.3. The process T = {Tt,t > 0) is a non negative Q-martingale. 

This fact allows to define the process := {Zj^ ,t > 0) absolutely continuous with respect to 
Z := {Zt,t > 0) on each Qt, t >0, with Radon Nykodim derivative: 

yAeGt, nz" eA) = E{lA{Z)Tt). 

We deduce from Lemma 13.21 and Lemma 13.31 the following Theorem. 

Theorem 3.4. We have that: 

(a) The limit of the empirical measure: 

1 ^ 

n=l 

exists a.s. 

(b) The process (Zf := Y/'R^,t > 0) %s distributed as {Z^,t> 0). 

Remark 3.5. We may interpret Theorem 13.41 as follows. The effect of the additive /i-transform 
factorizes in two parts, according to the decomposition of the Radon Nikodym derivative: 

Zt{ht) = Yt RM. 

The first term Yf induces a size bias of the total population size Z^{1) = Y^, see formula f[TH|) . 
whereas the second term Rt{ht) forces the first level particle to follow an /i-transform of P, see 
formula (fT9l) . 

The sequence {X^{n), G N) is not exchangeable in general, which contrasts with the initial 
sequence (Xt(n),n G N). The following Proposition shows that, loosely speaking, the first level 
particle is precursory. 

Proposition 3.6. Conditionally on {R^ = fi}, Xi{t) is distributed according to: 

F{X^il) G dx) = ^/i(rfx), 

and {X^{n))n>2 is an independent exchangeable random sequence with de Finetti's measure fi. 
3.2.2. Proofs. 

Proof of Lemma \3.2[ It is enough to observe that, by construction, the law of {Y^,U'^,X^) 
is absolutely continuous with respect to the law of (F, U, X) on J^t, with Radon Nykodim 
derivative given by: 

YpRoiho) ho{X^{l)) Yt m{0) m{t) 
MYoRoiho)) i?o(M Yom{t)ho{X^{l))m{0) 



vA G J-,, p((F^ u\x'^) eA) = E{ jrzj,. V'. . ^ lui i Hi uy, u,x) 



This also yields (the obvious fact) that {St,t > 0) is a J-"- martingale. □ 
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Proof of Lemma \3.^ From Lemma \2A\ since Qt C J^t and S* is a J-'-martingale, the projection 
E(S'fl^t) is a ^-martingale. We also have: 

where we used that Xtil) has law Rt conditionally on Qt for the third equality. Thus (T^, t > 0) 
is a ^-martingale. □ 



Proof of Theorem 3.4 ■ From Lemma [3.2[ the law of X'^ is absolutely continuous with respect 
to the law of X. The existence of the a.s. limit of the empirical measure of X'^ follows from 
that of X (but not the exchangeability of the sequence) and yields point (a). We prove point 
(b) now. Take A & Qt- 

¥{Z'' eA)=¥.{StlA{Z)) 

= ^{¥.{St\Gt) 1a{Z)) 
= P(Ti 1a{Z)) 
= ¥{Z" G A), 

where we use Lemma 13.31 at the third equality and the definition of for the last equality. 

□ 

Proof of Proposition \3M Let n G N be fixed, and let (0i)(i<i<n) be a collection of bounded and 
measurable functions on E. 



2<i<n 



E(Z„(A„)) 'Vfl,(fc,), 



2<i<n 



where we use Lemma [3.21 at the first equality, the de Finetti Theorem at the third equality, and 
Theorem 13.41 at the last equality. Since functions of the type Y[i<i<n^i characterize the law of 
n-uple, this proves the Proposition. □ 
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3.3. Applications. Overbeck investigated in [13] /i-transform of measure valued diffusions, 
among which the Dawson Watanabe process (with quadratic branching mechanism) and the 
Fleming Viot process (which is the GFV process for z/ = 0) using a martingale problem ap- 
proach. He also provided a pathwise constructions in the first case, see [13]. We shall see in this 
last Section how Theorem 13.41 applies in both cases and sheds new light on Overbeck's results. 

3.3.1. Generalized Fleming Viot processes. Recall the Generalized Fleming Viot process (with 
mutation) is the process Z constructed in Section [XT] when setting: 

. y = 1, 

• [/ is a subordinator with jumps no greater than 1, 
It also corresponds to the process R introduced in Section [5TT] when allowing for mutations. We 
denote by the Laplace exponent of the subordinator U : 

0(A) = c\+ [ (1 - e-^")z/^(c/x) 

J {0,1] 

where c > and the Levy measure u'^ satisfies J^q^^x (dx) < oo. The genealogy of the 
lookdown particle system is by construction described by the A-coalescent of Pitman [T3] . The 
finite measure A is related to u and c by x~^A|(o,i]((ix) = v{dx) and A{0} = c, and may be 
recovered from the characteristics of through the identity: 

/ g{x)K{dx) = cg{Q) + / (7(v^)z/^(rfx), 

J[0,1] >/(0,l] 

see Section 3.1.4 of [3]. 

Since li = 1, F is a martingale and we may apply results of Section 13.2.11 for any non- 
negative space time harmonic function {ht{x),t > 0,x & E) for the spatial motion P, that 
is any function such that {ht{^t),t > 0) is a non-negative martingale where ^ stands for the 
canonical process under P. Notice the construction of the particle system simplifies here 

since (f/^F'') ^'"= {U,Y) = ([/,!). 

Overbeck suggested in [T3] that in the particular case of the FV process, an additive h- 
transform looks like a FV process where ''the gene type of at least one family mutates as an 
h-transform of the one particle motion". This suggestion was made ''plausible" by similar 
results known for superprocesses, see [13] or the next Subsection, and a well known connection 
between superprocesses and Fleming- Viot processes which goes back to Shiga [20] • We did not 
attempt to derive the additive /i-transform for GFV processes in this way, since the connection 
between superprocesses and GFV processes is restricted to stable superprocesses and Beta GFV 
processes, see Birkner et. al [2]. 

Theorem 13 .41 allows us to see at first glance that the family which "mutates as an h-transform" 
is the family generated by the first level particle in the lookdown process. Other advantages 
of Theorem 13.41 over the martingale problem approach is that it provides a pathwise approach 
and also applies for GFV processes. 

Remark 3.7. If z/^ = (or, equivalently, A{dx) = A{0} 6o), the truncated processes obtained 
by considering the first particles: 

Z^{dx) Yl ^Xt{:n){dx) and Z^'^{dx) := ^ J] ^x^{n){dx) 

l<n<N l<n<N 
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correspond respectively to the Moran model with particles (see [3]) and its additive h- 
transform. This proves our approach is robust, in the sense that we can also consider discrete 
population. 

Finally, we may interpret the /i-transform as a conditioned process. For fixed s > 0, the 
additive /i-transform of the GFV process on [0, s] may be obtained by conditioning a random 
particle chosen at time t, t large, to move as an h-transform. For other conditionings on 
boundary statistics in the context of measure valued branching processes, we refer to Salisbury 
and Sezer [T8] . 

3.3.2. The Dawson Watanabe superprocess. Recall a continuous state branching process is a 
strong Markov process characterized by a branching mechanism ip taking the form 

(21) ^{\) = ]-a^\^ + 13\+ f (e-^"-l + AMl„<i)z/^(rfM), 

^ J(0,oo) 

for a Levy measure such that /(gj^-)(l A u'^)v^ {du) < oo, /3 G M, and G We will 
denote it CB(?/') for short. More precisely, the CB(?/') process is the strong Markov process Y 
with Laplace transform given by: 

E(e-^^* 1^0 = x) =e-^"(^'*\ 

where u is the unique non-negative solution of the integral equation, holding for all t > 0, 
A > 0: 

(22) M(A,t)+ f ds ^ {u{\ s)) = \. 

Jo 

We assume that ip'{Q+) > — oo, so that the CB{iIj) has integrable marginals, and 

(rie'^'(°+)*,t>0) 

is a martingale. The Dawson Watanabe process with general branching mechanism given by 
is the measure valued process {Zt,t > 0) constructed in Section \3A] when: 

• {Yt,t> 0) is a CB. 

• {Ut,t > 0) is the quadratic variation process of Y, Ut = [Y](t). Therefore, AUt = 
{AY,f < Y,\ 

Since the process (F^ e'^'*^^"'"^*, t > 0) is a martingale, we may apply our results for any non- 
negative function {ht{x),t > 0, x E E) such that {ht{^t) e""^ > 0) is a martingale. 
We now link our results with the literature: 

(1) When y is a subcritical CB process, meaning that ^'{0+) > 0, setting m{t) = e"'^''-*^"''-'* 
and ht{x) = e'^ '-^"'"-'*, we recover from Theorem 13.41 part of the Roelly & Rouault [I6] 
and Evans [3] decomposition. This /i-transform may be interpreted in that case as the 
process conditioned on non extinction. 

(2) When F is a critical Feller diffusion and P the law of a Brownian motion, if we assume 
we are given ht{x) a space time harmonic function for P and set m(t) = 1, we get from 
Theorem 13.41 the decomposition of the /i-transform of the Dawson Watanabe process 
provided by Overbeck in |13j . 
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We identified in Remark 13.51 two effects of tlie additive /i-transform: tlie total population is 
size biased and the first level particle follows an /i-transform of P. We shall concentrate on the 
first effect, and explain how a "spinal" decomposition may be partly recovered from Theorem 
13.41 Lemma 13.81 identifies the size biased total mass process = Z^{1) with a branching 
process with immigration, and Lemma [3.91 recognizes the first level particle as the source of the 
immigration. 

Let be the Laplace exponent of a subordinator. Recall a continuous state branching process 
with immigration with branching mechanism ip and immigration mechanism 0, CBI('0,0) for 
short, is a strong Markov process {Y^,t > 0) characterized by the Laplace transform: 

£^g-Ay/ |y-j _ ^-j _ ^-xu[X,t)-J^^ds 0(«(A,s)) ^ 

We recall for the ease of reference the following well known lemma, and we stress that this 
Lemma also holds in the supercritical case ip'{0+) < 0. 

Lemma 3.8. The process defined by with m(t) = e"^'*^*^"*"^* is a CBI(il},(j)) with immi- 
gration mechanism given by 0(A) := — ^/''(0+). 

The proof is classical and relies on computation of the Laplace transforms. Notice that in the 
case where the CB process Y extincts almost surely, the CBI process Y^ may also be interpreted 
as the CB process Y conditioned on non extinction in remote time, see Lambert [TU] . 

The total mass process Y^ = Z'^{1) is from Lemma 13.81 a CBI process whereas Y is by 
assumption a CB process. We may thus wonder "who" are the immigrants in the population 
represented by the particle system X^. The following Lemma shows that the offsprings of the 
first level particle are the immigrants. Recall ji refers to the first level sampled in the lookdown 
construction. Let us denote Ji(s) instead of ji for indicating the dependence in s. 

Lemma 3.9. The process (X]o<s<t ^^1^ ■'-{ii(s)=i}' ^ — O) ^■^ ^ pure jump subordinator with Levy 
measure up^ {du). 

Proof. By assumption, the process F is a CB{%1)) and from Lemma [3^ Y^ is a CB\{il), 0). From 
the Poissonian construction of CBI, we have that the point measure 

X] ^{s,AYi^)ids,du) 

0<s<t 

has for predictable compensator 

ds {Y^tu^{du) + uu^{du)). 

The expression of the compensator may be explained as follows. The term ds Y!}_v'^ [du) 
comes from the time change of the underlying spectrally positive Levy process, called the 
Lamperti time change (for CBs). The term ds uv^ [du) is independent of the current state of 
the population and corresponds to the immigration term. Then, conditionally on the value of 
the jump dsY^ = u, the event {ji{s) = 1} has probability 

u u 

Y^ ~ Y^_ + u 

independently for each jump. Therefore, the predictable compensator of the point measure 

^(s,Ayh)(c?s, (iM)l{jj(5)=i} 

0<s<t 
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IS 



u 



ds ( — r (^s-^ {du) + uu idu)) = ds uu (du), 

This ends up the proof. □ 

Remark 3.10. Understanding the action of the continuous part of the subordinator requires to 
work with the discrete particle system generated by the first TV particles. Namely, it is possible 
to prove that the family of processes 



H m<^<N.3^{ s)<N} ^ 

l{iiW=ij2(s)<Ar},* > U 



N 

yO<s<t 

converges almost surely as — ?■ oo in the Skorohod topology towards a subordinator with 
Laplace exponent </>. 
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